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ABSTRACT 

We present Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) images of the 
merger remnant NGC 7252. In particular, we focus on the surface brightness profiles 
and effective radii R c g of 36 young massive clusters (YMCs) within the galaxy. All the 
clusters have masses exceeding 10 5 M and are, despite the 64 Mpc distance to the 
galaxy, (partly) resolved on the HST images. Effective radii can be measured down to 
~ 2.5 pc, and the largest clusters have R e s approaching 20 pc. The median R e s of our 
sample clusters is ~ 6 — 7 pc, which is larger than typical radii of YMCs (~ 2.5 pc). 
This could be due to our sample selection (only selecting resolved sources) or to an 
intrinsic mass-radius relation within the cluster population. We find at least three 
clusters that have power-law profiles of the Elson, Fall, & Freeman (1987, "EFF") type 
extending out to > 150 pc. Among them are the two most massive clusters, W3 and 
W30, which have profiles that extend to at least 500 and 250 pc, respectively. Despite 
their extended profiles, the effective radii of the three clusters are 17.2, 12.6 and 9.1 pc 
for W3, W26 and W30, respectively. We compare these extended profiles with those 
of YMCs in the LMC (R136 in 30 Dor), the Antennae galaxies (Knot S) and in the 
nearby spiral galaxy NGC 6946. Extended profiles seem to be a somewhat common 
feature, even though many nearby YMCs show distinct truncations. A continuous 
distribution between these two extremes, i.e. truncated or extremely extended, is the 
most likely interpretation. We suggest that the presence or absence of an extended 
envelope in very young clusters may be due to the gas distribution of the proto-cluster 
giant molecular cloud, in particular if the proto-cluster core becomes distinct from the 
surrounding gas before star formation begins. 

Key words: galaxies: individual: NGC 7252 - galaxies: star clusters: general. 



1 INTRODUCTION 

NGC 7252 is a prototypical merger remnant that hosts one 
of the largest populations of young massive star clusters 
known (Whitmore et al. 1993; Miller et al. 1997, hereinafter 
M97), including two clusters with stellar masses exceeding 
1O 7 M (Maraston et al. 2004, hereinafter M04; Bastian 
et al. 2006, hereinafter B06). While there is ongoing star- 
formation within a central disk of gas (Schweizer 1982) the 
cluster population as a whole is dominated by an extended 
halo population of YMCs with an age of ~ 400 ± 100 Myr 
(e.g., Schweizer & Seitzer 1998, hereinafter S98). 

Studies of luminosity profiles of YMCs have shown that 



they often differ from those of classic globular clusters. The 
latter are generally well fit by King (1962) profiles that dis- 
play a truncation at large radii. This truncation is thought to 
be related to the tidal radius, inside which stars are bound 
to the cluster, while outside the stars move according to 
the galactic potential. Many YMCs, on the other hand, do 
not display such a truncation, but rather are well fit by ex- 
tended power-law envelopes. This was quantified by Elson, 
Fall, & Freeman (1987, hereinafter EFF) for YMCs in the 
Large Magellanic Cloud. These authors fit profiles of the 
form I(r) = Jo(l + 7 ,2 /a 2 )~ 7/ ' 2 , where r is the distance from 
the cluster centre and a is a characteristic radius. For r ^> a 
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the EFF profile becomes a simple power-law. EFF profiles 
have been shown to also provide good fits to YMCs in M33 
(San Roman et al. 2012) and in spiral galaxies more gener- 
ally (Larsen 2004). 

Observations of young cluster systems have shown that 
the clusters themselves are often grouped into larger struc- 
tures, cluster complexes, with radii of tens to hundreds of 
parsecs (e.g., Zhang et al. 2001; Larsen 2004). These com- 
plexes often appear to be centrally concentrated, hosting 
a dominant cluster in the centre, with a similar EFF-type 
power-law decrease in surface brightness with increasing ra- 
dius (Bastian et al. 2005). The outer envelope of the com- 
plexes is, at least partly, made up of smaller individual clus- 
ters. Fellhauer & Kroupa (2005) carried out iV-body model 
simulations of a system of massive star clusters distributed 
as a Plummer sphere, assuming that the clusters are in virial 
equilibrium, and found that the clusters merge within a few 
hundred Myr. The remnant formed in their simulations has 
a massive central cluster and an extremely extended outer 
halo, leading the authors to suggest that the massive cluster 
W3 in NGC 7252, may have formed in such a manner. 

Interested by the extreme nature of the YMCs in 
NGC 7252, i.e. their high masses and relatively large ages 
(~ 400 Myr), we have obtained HST/WFC3 images of 
NGC 7252 to study the luminosity profiles of these clusters 
in detail. Despite the relatively large (64 Mpc for Ho = 75 
kms -1 Mpc -1 ) distance to the host galaxy, a large number 
of the clusters are (semi) resolved. We were able to trace the 
profiles of three of the most massive clusters out to distances 
of more than 150 pc from the cluster centres. Specifically, 
our main goal was to see whether the cluster profiles are 
truncated and, if not, how far they extend. We were also 
interested in studying whether the presence or absence of 
any extended envelope might correlate with the projected 
distance of clusters from the remnant's centre, which would 
be expected if the tidal truncation of clusters in the galaxy 
potential would already have had an observable impact. 

The present paper is organized as follows. In § [2] we 
present the data and methods used to create the point 
spread function, and then outline the steps taken to gen- 
erate the cluster sample. In §[3] we introduce the techniques 
used to measure the effective radii of the clusters and study 
the luminosity profiles of three massive clusters (W3, W26 
and W30) in detail. In §[4] we study the distribution of R c g 
for the entire population and investigate the possibility of 
a mass-radius relation for massive clusters. In §[5] we inves- 
tigate the relation between the effective radius and galac- 
tocentric distance of the clusters in our sample. In § [6] we 
compare the extended profiles of clusters in NGC 7252 with 
those of other known clusters and discuss the possible origin 
of such extended envelopes. Finally, §[7] presents our conclu- 
sions. 



2 OBSERVATIONS 
2.1 WFC3 observations 

Broad-band images of NGC 7252 were obtained with 
HST/WFC3 on 2010 July 29 through the F475W (SDSS 
g) and F775W (SDSS i) filters (GO-11554, PI: N. Bastian). 
The filters were chosen to provide coverage in the blue and 



red part of the optical, and to be wide in order to maximise 
the throughput and signal-to-noise ratio of the brightness 
profiles. 

The total exposure time through each filter was 2120 s, 
split into four 500 s exposures plus two short, 60 s exposures 
to avoid saturation of the bright central regions of each clus- 
ter and the galaxy. The individual exposures in each filter 
were taken with small spatial offsets ("dithers") between 
one another to allow improved spatial sampling of the point 
spread function (PSF) as well as the elimination of hot pix- 
els and cosmic-ray hits during image combination. Prior to 
running the PyRAF/stsdas task MultiDrizzle to produce 
a combined image for each filter, we created sky variance 
maps for each individual exposure for the purpose of de- 
riving weight maps for the final image combination. These 
maps were constructed from the sky values of the individ- 
ual exposures, the flatfield reference file, the dark current 
reference file, and the read-out noise values as listed in the 
image headers. As such, these sky variance maps contain all 
noise components of the image except for the Poisson noise 
associated with the sources on the image. 

The final run of MultiDrizzle was performed by shrink- 
ing the input pixels by 20% at the stage where input pix- 
els are drizzled onto the output image grid (i.e., a value of 
0.8 was adopted for the f inal_pixf rac parameter), and we 
chose an output image pixel size of O'.'028/pixel. These var- 
ious parameters were selected after extensive experimenta- 
tion, and they allow a good match to the degree of subsam- 
pling induced by the dither pattern we used. Saturated pix- 
els in the 500 s exposures were replaced by the corresponding 
pixels in the short exposures while running MultiDrizzle. 
This was achieved by setting the appropriate data quality 
flag for the affected pixels of the long exposures. 

Additionally, we used WFC3 F336W (U) images (GO- 
11691, PI: P. Goudfrooij) that will be presented in more de- 
tail in Goudfrooij et al. (in prep.). In brief, the total F336W 
exposure was 3050 s, and the images were reduced in much 
the same way as those described above. These U images 
were not used for size or profile determination, but only for 
photometric selection. 



2.2 Other imaging 

In addition to the WFC3 images of NGC 7252, we also used 
archival data to study the luminosity profiles of other YMCs. 
In particular, we used VXT/HAWK-I near-infrared .Ks-band 
images of 30 Doradus to study the outer profile of the central 
cluster, R136 (programme ID 60.A-9283). The data were 
taken as part of the commissioning of the instrument and 
are presented in more detail in Campbell et al. (2010). 

For the YMC in NGC 6946, we used the #ST/WFPC2 
1^-band images presented in detail by Larsen et al. (2001). 
Finally, for the study of Knot S in the Antennae galaxies 
we used HST/ ACS l/-band images taken from the Hubble 
Legacy Archive. These images are presented in more detail 
in Bastian et al. (2009) and Whitmore et al. (2010). 

2.3 Creation of a model PSF 

In order to get accurate measurements for the profiles and 
radii of the target clusters we required an accurate model of 
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the point-spread function (PSF) of our data. We did this in 
the following two ways. 

The first was to select sources of various brightnesses 
throughout the image that appeared to be unresolved (i.e., 
they all had the same FWHM) and use them in the IRAF 
task PSF to generate an empirical PSF. For this we used 14 
sources located throughout the image. The disadvantage of 
this method is that it will blur any effect of the PSF varying 
across the chips. 

A second set of PSFs was created from the grid of em- 
pirical WFC3 PSFs assembled by Jay Anderson (hereafter, 
ePSFs; see Anderson & King 2006 for the ACS version). 
Since the latter PSFs are created for individual *_flt.fits 
images rather than for images combined with MultiDrizzle, 
we wrote a package based on Maurizio Paolillo's Multik- 
jijEI suite of scripts (modified to work for WFC3 images). 
The package creates empirical PSFs at the location of each 
GC candidate in blank versions of each individual input 
_f It. fits WFC3 image, and then combines those files by 
MultiDrizzle in the exact same way as the final NGC 7252 
images. These PSFs were then subsampled by a factor 10 
in order for them to be used appropriately within ISHAPE 
(see below), using polynomial interpolation. 

We compared the empirically derived PSFs with the 
PSFs in the ePSF grid by deriving the radii of the spec- 
troscopically confirmed clusters (see § 12.41) using ISHAPE 
(Larsen 1999; see also § 13. II below) with both sets of PSFs. 
We did not find any systematic differences between the re- 
sults using either of the PSFs, with measured FWHMs being 
the same to within ~ 10%. This difference is smaller than 
that due to the choice of the radius over which to measure 
the size, or the choice of the type of profile to fit (King vs 
EFF). For the results presented below we adopt the ePSF 
grid. 

For each of the methods we treated the two filters 
(F475W and F775W) independently. 

2.4 The cluster sample 

The main goal of this work was to measure the surface- 
brightness profiles (SBPs) and R e g of the clusters within 
NGC 7252, focussing on the ~ 400 ± 100 Myr old popula- 
tion formed during the galaxy merger, i.e., clusters outside 
the central 12"-diameter star-forming disk (M97, S98). This 
population is known to include extremely massive clusters 
(M97, S98, M04, B06), even though the population is con- 
sistent with a power-law mass function of index —2 (B06). 

In order to study the SBPs of the clusters, high S/N 
is required. We selected ~ 100 source candidates (18 < 
F475W < 25 mag) outside the central disk (galactocentric 
radii > 4.2", ~ 1.3 kpc), extending out to 108" (~ 33 kpc) 
from the galaxy center. We excluded candidates with nearby 
neighbours that would complicate the SBP analysis. In order 
to remove background galaxies we measured the magnitudes 
and colours of each of these candidates from the F336W, 
F475W , and F775W images, using an aperture of 10 pixels 
in radius. The resulting colour-colour diagram is shown in 
Fig. [T] where the spectroscopically confirmed clusters (S98) 
are shown in red. The dashed box shows the colour region 

1 http: //www. na.infn.it/~paolillo/Software. html 



where we selected our sample. Additionally, we show the 
evolutionary track of a GALEV simple stellar population 
(SSP) model of solar metallicity as a solid line (Kotulla et 
al. 2009). The limiting ages of the model clusters that pass 
our colour selection are 250 Myr and 1 Gyr, as marked. Note, 
however, that the observed cluster colours are concentrated 
around 400-500 Myr. The applied colour selection removes 
background galaxies and effectively limits our study to the 
intermediate age (~ 400 Myr) cluster population. Using our 
colour selection criteria, we found 52 candidate clusters. 

We then fit these candidates with ISHAPE, adopting an 
EFF profile with 7 = 3.0, and measured their concentration 
index (CI; see § 13.11 below). We also performed aperture 
photometry for each of the candidates with aperture radii 
ranging from 1 to 30 pixels. Examples of the resulting growth 
curves are shown in Fig. [2] 

Using these two methods, we selected objects that ap- 
peared extended in each method. Figure shows a com- 
parison of the measured CIs and R e « for an EFF profile 
with 7 = 3.0. The dashed lines mark limits where we can 
resolve sources. Red triangles represent confirmed clusters 
within NGC 7252 based on optical spectroscopy and agree- 
ment with the systemic velocity (S98). Additionally, we per- 
formed resolution tests by adding artificial clusters (with 
known radii) to the images which were analysed in the same 
way as discussed above, in order to confirm the adopted 
limits. 

As Fig. [3] clearly shows, all of the confirmed clusters 
from S98 are resolved. Also, of the total sample of 52 can- 
didate clusters 36 are resolved in both the ISHAPE and CI 
analyses. However, we note again that our cluster sample is 
not complete in any way, as only bright and resolved clusters 
were selected. 

Table[2]gives the IDs, magnitudes, CIs, estimated effec- 
tive radii (adopting an EFF profile with 7 = 3.0) and galac- 
tocentric distances of all the clusters that pass our colour 
selection and appear resolved using both ISHAPE and the 
CI method. 



3 SURFACE-BRIGHTNESS PROFILES AND 
RADII OF THE MASSIVE CLUSTERS 

3.1 ISHAPE and the Concentration Index 

Once the cluster candidates were selected we used three 
methods to study their profiles and effective radii. 

First, for the brightest clusters (e.g., W3, W26 and 
W30) we performed a detailed analysis of their SBPs with 
the intention to constrain the profile shape, effective radius 
and extent to which their envelopes reach. These will be 
discussed individually below. 

As a second method, we used the ISHAPE algorithm 
(Larsen 1999) with a fixed profile shape (an EFF profile 
with 7 = 3) to measure each cluster's FWHM, which in 
turn can be used to derive its R B a. The fitting radius was 
chosen to be 20 pixels (= 0.56" ~ 175 pc), although we did 
not find any systematic effects as a function of the fitting 
radii. 

As a third method, we measured for each cluster the 
concentration index CI (e.g., M97; Whitmore et al. 2010), 
which is defined as the magnitude difference between mea- 
surements with a 1 pixel and a 3 pixel radius aperture. We 
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Figure 1. The colour-colour diagram of isolated sources in 
NGC 7252. Red points represent clusters that have been spectro- 
scopically confirmed, while the dashed box indicates the colour 
selection applied. All sources within the box are considered clus- 
ter candidates. The solid line represents the GALEV SSP models 
for solar metallicity, and the limiting ages of the models that pass 
our colour criteria are given in the panel. 



Figure 3. Estimated effective radii R e g determined by using 
ISHAPE (adopting an EFF profile with 7 = 3) are plotted ver- 
sus concentration indices CI(1— 3) measured with apertures of 
one and three pixel radius. The vertical dashed line marks the 
FWHM of 0.2 pixels used for the adopted profile, while the hor- 
izontal dashed line denotes the cut applied between point-like 
and extended sources. The spectroscopically confirmed clusters 
are shown as (red) triangles. Only sources that passed our colour 
selection arc shown. 




1 10 

Aperture (pixels) 



Figure 2. Curves of growth, normalised to an aperture of 1 pixel, 
as a function of aperture radius for a sample of sources in the 
WFC3 FOV (before any colour selection has been applied). The 
inset shows a zoom in on the end of the distribution, with con- 
firmed clusters and likely stars labelled. Confirmed clusters are 
shown in black with filled circles, while other sources in the field 
are shown as dashed grey lines. 



also experimented with different aperture combinations and 
found consistent results. 

In all plots of the SBPs to be shown below, we have 
added a constant to the derived SBP so that the surface 
brightness becomes mag arcsec -2 in the first bin. When 
profiles of more resolved clusters are shown for comparison, 
we shift these to have a value of mag arcsec -2 at the same 
physical radius as that for the NGC 7252 clusters. 



3.2 W3 

3.2.1 Profile and effective radius 

Cluster W3 is an enormously massive cluster of ~ 8xl0 7 M 
(S98, M04). In order to determine its profile and effective ra- 
dius we fit various profile types with ISHAPE, letting the 
FWHM as well as the index 7 (for an EFF profile) or con- 
centration parameter c (for a King profile) vary. We tested 
multiple initial guesses (index and concentration) for the al- 
gorithm and found good convergence. 

Once the best fitting profiles were found for each profile 
type, we created artificial clusters of that profile type and 
radius, with the same brightness as W3 (within 40 pixels) 
and added them to the images (7 artificial clusters for each 
profile type) at the same galactocentric distance as W3, but 
at different position angles. We then directly compared the 
SBPs of the models with W3. Example profile comparisons 
are presented in Figs. [4] and [5] for King and EFF profiles, 
respectively. 

As can be seen from Fig. [4] the profile of W3 does not 
show any signs of truncation out to where the background 
was taken at 2.2" (« 700 pc). Comparing with the best fit- 
ting King model profiles of c = 100 (top panel) and c = 300 
(bottom panel), it is clear that — if present-such a trunca- 
tion would have been observed. If we allow c to vary, values 
of ~ 500 are preferred. 

In contrast, EFF model profiles, which do not have any 
truncation, provide excellent fits to the data. Figure [5] shows 
two examples: the lowest-;^ 2 profile (within a fitting radius 
of 20 pixels = 0.56" ~ 175 pc) is shown in the top panel, 
while a similar profile with 7 = 2.6 is shown in the bottom 
panel. Both model profiles reproduce the data quite well, 
and we adopt a best-fit index of 7 = 2.4 ± 0.2 for the full 
cluster profile. 

The corresponding effective radius is R c g — 17.2l 2 P c 
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Figure 4. The surface brightness profile of W3 (red circles) along 
with the best fitting King profile (green triangles) with c = 100 
(top) and c = 300 (bottom). The fits were carried out within 
0.56" and the model profiles were measured on artificial clusters 
added to the image at the same galactocentric distance as W3. 
The dashed vertical line marks a 0.7" radius (~ 215 pc) for ref- 
erence. Note that a truncation in the profile in either case would 
be observable. If the underlying profile is best described by a 
King profile, then in order for the truncation to not be seen the 
concentration parameter would need to be c > 500. 



for both the F475W and F775W images. The reason for 
the larger error on the high end is that for shallower indices 
(approaching 7 = 2) R e s becomes undefined. Hence, values 
of 7 near 2 lead to large uncertainties in the determination 

Of iZeff- 



3.2.2 Dynamical mass 

Using the best-fit model profile (EFF with 7 = 2.4) we 
estimate an effective radius of Y1.2 J [\ pc for W3. This 
value is very similar to that found by M04, namely R e g = 
17.5 ± 1.8 pc, based on lower-resolution WFPC2 observa- 
tions. Hence, we confirm the dynamical mass estimate of 
8 ± 2 x 10 7 M by M04, making W3 the most massive star 
cluster presently known. 



Figure 5. Similar to Fig. [4] but now the model profiles are for 
clusters with an EFF profile. The top panel shows the best fitting 
(lowest x 2 ) profile, while the bottom panel shows a profile with 
7 = 2.6. Such power-law like profiles provide a better fit to the 
outer envelope of W3 than King profiles. 



3.3 W30 

The next most massive cluster in the NGC 7252 popula- 
tion is W30, with a measured dynamical mass of 1.6 ± 0.3 x 
1O 7 M (B06). W30 is ~ 1.6 mag fainter in the V-band 
than W3, resulting in a lower S/N, especially in the outer re- 
gions. Additionally, the analysis is complicated somewhat by 
the presence of two neighbouring sources at 0'.'92 and l'.'04. 
However, these sources are significantly fainter than W30 on 
both the F475W and F775W images and can, therefore, be 
masked out by the ISHAPE software during the fitting. 

The measured SBP of W30 is shown in Fig. [6] as solid 
(green) circles. One can readily see by eye that W30 is more 
compact than W3, although it, too, features an extended 
profile. This profile can be traced out to at least 250 pc 
(0.8") from the cluster centre. 

Carrying out an analysis similar to that described above 
for W3, we estimate values of R c g — 7.45 ± 0.22 pc and 
10.8 ± 1.4 pc for W30 on the F475W and F775W im- 
ages, respectively. The straight average of these two values 
is R c g — 9.1 ± 1.7 pc, in excellent agreement with the value 
of 9.3 ± 1.7 pc found by B06. 
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Figure 6. Intcrcomparison of the profiles of three extended clus- 
ters, shown along with the profile of a foreground star in the same 
field of view. Note that W3 and W30 extend to at least 500 pc 
and 250 pc, respectively. The fact that the foreground star used 
as a comparison crosses the extended profiles of W3 and W30 is 
due to its faintness, making it difficult to reliably trace its profile 
past 8 magnitudes below its central surface brightness. 



Table 1. Cluster ID, effective radius, and best fit profile index 
(assuming an EFF profile) for the clusters discussed in §[3] 



ID 


R-cff 


7 (EFF profile) 


Extent (pc) 


W3 

W30 

W26 


17.2 +b 2 

Q 1+2-3 
a,1 -2.3 

12.611;? 


2 4+ ' 2 

2 4+ ' 2 
z -^-0.2 
n 7+0.2 
z - '-0.2 


> 500 

> 250 
150 



The index of the best-fit EFF profile for W30 appears 
to be similar to that for W3, with 7 = 2.50 ± 0.04 and 
7 = 2.36 ± 0.06 on the F475W and F775W images, respec- 
tively. There is some indication that W30 is slightly elon- 
gated (minor/major axis = 0.84) on the blue image, but no 
trace of this is found on the red image. 

3.4 W26 

Another cluster that stands out in Fig. [3] is W26, the profile 
of which is shown in Fig. [5] This cluster is ~ 0.9 magnitudes 
fainter than W30 (M97), making it more difficult to measure 
the outer parts of its profile. This may perhaps explain the 
relatively abrupt truncation observed at 0'.'5 (~ 150 pc), but 
alternatively the truncation could also be real (§ 13. 6[) . 

We measure R c a — 11.4 pc and 13.85 pc in the blue and 
red bands, respectively, which yields a straight average of 
R cB = 12.6 ± 1.7 pc. The EFF-profile index of W26 appears 
to be steeper than that of W3 or W30, with 7 = 2.9 and 2.5 
in the blue and red bands, respectively. No evidence for any 
elongation is seen for W26 in either band. 

3.5 Other massive clusters 

In addition to the three clusters with extended profiles dis- 
cussed above, we attempted to derive the profile shapes for 
a number of other bright massive clusters in NGC 7252. Due 



to the lower S/N of these clusters we did not create artifi- 
cial cluster models and add them to the image to investigate 
their profiles directly; instead, we simply used ISHAPE to 
estimate 7, adopting an EFF profile. In order to test the 
reliability of the estimated parameters, we ran the fits with 
three initial guesses for 7, namely 7 = 2.0, 4.0 and 6.0, and 
required the estimated values to converge to the same value 
for all initial guesses. 

In this way, we were able to estimate 7 values of 2.8, 2.6, 
3.0 and 3.0 for W6, S105, S114 and W10, respectively. These 
values are near to, although slightly steeper than, those de- 
rived for W3, W30 and W26. Due to the similarity between 
the estimated 7 values and that of 7 = 3.0 adopted for es- 
timating 7? c ff for the full sample, the radii of these clusters 
do not differ significantly from those given in Table [2] 



3.6 Tidal radii 

Given the extended envelopes around W3, W30 and W26, it 
is insightful to relate their extent to an estimate of the tidal 
radius for each cluster. The tidal radius can be estimated 
via 

where M c i is the mass of the cluster, Rq the true distance 
from the galaxy centre, and Vg the circular velocity within 
the galaxy at that distance (von Hoerner 1957). This esti- 
mate assumes that the clusters are on circular orbits, and we 
approximate the true distance with the projected distance. 
If the clusters are located further away from the galaxy cen- 
tre their tidal radii would be larger, and if the clusters move 
on radial or elliptical orbits (i.e., passing near the centre of 
the galaxy) their tidal radii would be significantly smaller 
than estimated via the above equation and approximations. 
Additionally, the above approximation may not be accurate 
if the three extended clusters are located at a true galac- 
tocentric distance near their projected one. In such a case, 
each cluster's extent would be a significant portion of the 
galactocentric distance, violating the initial assumptions of 
the approximation. 

For the masses of W3, W30 and W26 we use 8 x 10 7 M 
(M04), 1.6 x 10 7 M Q (B06) and 6 x 10 6 M (assuming the 
same age as W3 and W30), respectively. Additionally, we 
assume a circular velocity of ~ 200 kms -1 . With these val- 
ues we estimate the tidal radii of W3, W30 and W26 to 
be approximately 450 pc, 200 pc and 150 pc, respectively. 
Given the various uncertainties associated with these esti- 
mates, the observed extents of the profiles of W3, W30 and 
W26 given in Table 1 appear largely consistent with the 
estimated tidal radii. Specifically, the observed extents of 
the envelopes of W3 and W30 are at least equal to the es- 
timated tidal radii and may exceed them, suggesting that 
either the true galactocentric distances of these two clus- 
ters exceed the projected distances significantly or there has 
been insufficient time for significant tidal erosion. The near 
agreement between the measured extent of W26 and the es- 
timated tidal radius, on the other hand, suggests that the 
observed truncation may be tidal and real, rather than due 
to a measuring problem. Our main conclusion is that the ob- 
served large extents and/or tidal radii of these clusters are 
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Table 2. Magnitudes, concentration indices CI, effective radii, and projected galactocentric distances for 36 young massive clusters of 
NGC 7252 discussed in § [3] Note that all structural parameters were measured assuming an EFF profile with 7 = 3.0, which is why 
some of the values differ from those in Table [T] The observed magnitudes have not been corrected for foreground or internal extinction. 
The IDs beginning with 'W refer to Whitmore et al. (1993), 'S' refers to S98 and 'Mi' to M97. For clusters not previously catalogued 
we provide the coordinates in the final two columns. 



ID F336W F475W F775W CI R cff (F475W) R cS (F775W) Galactocentric RA and DEC 

(mag) (mag) (mag) (mag) (pc) (pc) distance (kpc) (J2000) 
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due, at least in part, to the high cluster masses. Hence, ex- 
tremely extended profiles may be a common feature of very 
high-mass clusters. 



Due to the large extent of W3, W30 and W26, we have 
also searched for tidal debris around each cluster. Figure [7] 
shows surface brightness maps (mag/arcsec 2 ) in the F475W 
image, of each cluster. All three clusters appear circular and 
no clear tidal debris is seen. From these images we can see 
that a the weakly varying background is not affecting our 
results (e.g., spiral structure, dust lanes, etc) and that bin- 
ning the profile of the clusters radially (i.e., averaging az- 
imuthally) does not hide structure. Tidal features, while not 
visible in the surface brightness images, may still be present; 
however the relatively bright background limits the possibil- 
ity of detecting such structures with the present data. 



4 RADIUS DISTRIBUTION AND THE 
MASS-RADIUS RELATION 

While we cannot fit SBPs for the other cluster candidates 
in our sample due to their lower S/N and blending with 
the background galaxy light, we have attempted to measure 
R a g for each of the candidates. For this we only use sources 
that are found to be resolved with both ISHAPE and the 
concentration index technique (see Fig. [3}, i.e. 36 of the 
52 cluster candidates identified through their colours (see 

§ESJ- 

We find that both the mean and median values of _R c ff 
are ~ 6-7 pc, independent of the filter used. This is sig- 
nificantly larger than for typical YMCs (~ 2.5 pc; e.g., 
Larsen 2004; Portegies Zwart et al. 2010). There are three 
possible explanations for this difference. The first is that 
the difference is real and that our subsample of clusters in 
NGC 7252 is representative of the size distribution of the 
cluster population as a whole within this galaxy. The sec- 
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Figure 7. Surface brightness maps (F475W) of the three clus- 
ters with profiles that extend beyond 150 pc. Each panel is 4" on 
a side (1.2 kpc) and no smoothing has been applied. The (rela- 
tively smooth) gradient in the background stellar brightness (from 
NGC 7252) can be seen in each panel. In the inner region of the 
clusters, each pixel has been set to the minimum surface bright- 
ness shown in the colour bar (i.e. the inner sections are not re- 
solved in this representation). Note that all three clusters appear 
circular in projection and that no tidal features are apparent. 
However, due to the relatively bright background, any possible 
tidal features would be difficult to discern. 




6.0 6.5 7.0 

log (Mass/Msun) 

Figure 8. The relation between mass and effective radius for the 
young massive clusters of NGC 7252. Masses were determined by 
assuming that all clusters in the sample have the same age (due 
to our colour selection), suffer no extinction, and can simply be 
scaled from their FA75W magnitudes. The (red) triangles mark 
spectroscopically confirmed clusters, while the open circles denote 
clusters that were used in the regression fit. The best-fit relation is 
shown as the solid line, while the dashed line shows the relation for 
luminous elliptical galaxies (Hasegan et al. 2005) for comparison, 
and the dash-dotted line marks the mean effective radius of 3.2 pc 
for old globular clusters. For masses above 1 x 10 6 M we find 
evidence for a mass-radius relation at the ~ 3a level. 



ond possibility is that we have biased our sample by only 
considering well resolved clusters. This is certainly affecting 
our sample at some level: for example, if we assign all the 
unresolved clusters an effective radius of 2 pc, the median 
of the distribution is < 5 pc. Finally, it has been suggested 
that there is an underlying mass-radius relation for clus- 
ters above 10 6 Mq (e.g., Hasegan et al. 2005; Kissler-Patig, 
Jordan, & Bastian 2006). Many of the clusters in our sample 
have masses above 10 6 Mq , so their large sizes may simply 
be a reflection of their high masses. 

In order to test this third possibility we searched for 
a mass-luminosity relation of the resolved cluster candi- 
dates that pass our colour selection. First, we translated 
each cluster's F475W magnitude to a mass by assuming 
that all clusters have the same age (400 Myr) — due to our 
colour selection — and solar metallicity (S98) and that they 
are not affected significantly by extinction. The result is 
shown in Fig. [S] The dashed line represents the extrapo- 
lated relation from luminous elliptical galaxies (see Hasegan 
et al. 2005), while the dash-dotted horizontal line marks 
the mean globular cluster radius, which is similar to our 
resolution limit. We performed a regression fit to the data 
above 1.0 x 10 6 M© (indicated by open circles) of the form: 
log(7? c ff /pc) = a + 6*log(M c i/M0), where M c \ is the cluster 
mass, and found a = -1.16±0.63 and b = 0.29T0.09; this fit 
is shown as a solid line. From this fit we find evidence for a 
mass-radius relation for clusters in NGC 7252 with masses 
above TO x 10 6 Mq at the ~ 3<r level. However, the exact 
significance of the correlation depends on the adopted limits 
(it is lower if a lower mass limit of 2.5 x 10 6 Mq is used). 
Additionally, if W3 is excluded from the fit, the significance 
of the mass-radius relation drops to ~ 1.5a. 



© 2012 RAS, MNRAS 000,rflfT2l 



Cluster profiles and sizes in NGC 7252 9 



5 THE EFFECTIVE RADIUS - 

GALACTOCENTRIC DISTANCE RELATION 

Figure [9] shows the measured effective radii of the resolved 
clusters in our sample plotted versus the projected galacto- 
centric distances Dgc- The radii were measured from the 
F475W image, using an EFF 7 = 3.0 model profile, and are 
given in Tabled Each cluster is represented by a small filled 
circle, while the large (red) points with error bars represent 
the mean values and standard deviations of R c a determined 
in six distance bins with six clusters each. The dashed line 
shows a linear regression fit to the individual clusters, which 
has a slope of 0.35 ± 0.20. 

In order to check the validity of the observed i? e ff - -Dgc 
relation we carried out a series of Monte Carlo simulations. 
First, we created a population of f 000 clusters with a true 
spatial distribution made to match the observed projected 
distribution of clusters in NGC 7252 (M97, esp. Fig. 19). 
Next we assigned a random _R c ff to each cluster, chosen from 
a gaussian distribution with a mean of 7 pc and a standard 
deviation of 4 pc. We then selected 36 clusters from the orig- 
inal 10 000 at random, matching their parameters to those 
of the observed population (1.85 kpc < -Dgc < 13.21 kpc, 
R e g > 2 pc). Finally, we carried out a linear regression on 
these clusters in the same way as was done for the obser- 
vations. This type of analysis implicitly assumes that the 
clusters move on circular orbits. 

One thousand realisations of the model were performed. 
Figure [TO] shows the resulting distribution of the slopes of 
linear regressions performed between R c a and galactocentric 
distance. The vertical dashed line marks the observed slope 
(shown in Fig. [SJ. We find that only ~ 5% of the random 
simulations resulted in a slope similar to or larger than that 
observed. 

We conclude that there is an intriguing relation between 
-Reg and the galactocentric distance of clusters in NGC 7252. 
However, the observed relation may be due to small-number 
statistics. Expanding the present study to other galactic 
merger remnants with massive clusters may allow a more 
definitive answer. Goudfrooij (2012) looked at the R e s - 
Dgc relation for massive clusters in the intermediate age 
(~ 3 Gyr - e.g., Goudfrooij et al. 2001) merger remnant 
NGC 1316. He found a similar relation as that reported 
here, that the mean effective radius increases as a function 
of projected galactocentric distance. This lends support to 
the notion that such a relation may be intrinsic to young-to- 
intermediate-age cluster populations in merger remnants. 
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Figure 9. The observed relation between R e g and projected 
galactocentric distance for the 36 resolved clusters from Table [2] 
The small filled circles represent the individual clusters, while the 
large (red) filled circles with error bars represent the mean R B ff 
and standard deviations of the clusters in distance bins contain- 
ing six clusters each. The dashed line shows the linear-regression 
fit to the unbinned data. 
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Figure 10. The distribution of measured slopes from 1000 Monte 
Carlo simulations of the NGC 7252 cluster population, assuming 
no intrinsic relation between R c g and galactocentric distance. The 
vertical dashed line marks the observed slope (see Fig- [9j ■ Only 
approximately 5% of the simulations yield regression slopes equal 
to or larger than the observed slope. 



6 DISCUSSION 

6.1 Comparison with Knot S in the Antennae 
Galaxies 

While W3 and W30 are two of the most massive star clusters 
known, they are not unique in having extended power-law 
envelopes. The Antennae galaxies host a number of massive 
clusters as well as extended "knots" of ongoing or recent 
star and cluster formation. The brightest such knot, known 
as 'Knot S' (Rubin et al. 1970; Whitmore et al. 1999), is 
located in the distorted disk of NGC 4038, away from the 
"overlap region" where the two galaxies' gas masses are cur- 



rently colliding; see Fig. 28 in Whitmore et al. (2010) for an 
JKT/ACS-HRC image of this knot. 

We used tfST/ACS-WFC V-band images (presented in 
detail in Bastian et al. 2009; Whitmore et al. 2010) to study 
the profile of Knot S (following Schweizer 2004). For the 
analysis we used the empirical PSFs described in Bastian 
et al. (2009). The profile of Knot S is shown in Fig. [TT1 as 
a solid (blue) line. We measured the profile of Knot S in 
the same way as for the NGC 7252 clusters discussed above, 
namely using ISHAPE, selecting an EFF model profile and 
fitting on the FWHM and index 7 of the profile. The initial 
estimate of Reg, for a fitting radius of 25 ACS pixels = 120 
pc, yields R c g — 18 pc, however with a shallow profile of 
7 ~ 2. Since the profile is so shallow, the effective radius 
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is not well constrained. If we fix a profile type and shape 
(EFF, 7 = 3.0) and fit both Knot S and W3 with it, we find 
that both fits result in i? c ff ~ 15-20 pc. Hence we conclude 
that both clusters have similar characteristic radii. Both also 
have envelopes of similar extent (> 450 pc). 

An interesting property of Knot S is its youth. Whit- 
more et al. (2010) estimate an age of ~ 6 Myr for the core of 
the knot, i.e. within 15 pc from the centre, through broad- 
band photometric techniques, and Whitmore et al. (1999) 
estimate a similar age, ~ 7 Myr, based on UV spectroscopy. 
Additionally, the significant amount of substructure (i.e. 
small subclusters) present in the halo of Knot S suggests 
that it is dynamically unevolved, i.e. the current distribu- 
tion is similar to the initial structure of the knot (although 
the inner regions may be dynamically more evolved). The 
subclusters in the envelope of Knot S show evidence for a 
range of ages from ~ 1 Myr to ~ 30 Myr. The oldest such 
subcluster is a massive cluster of ~ 1 x 1O 6 M0, seen as 
a bump in the SBP at R ~ 220 pc. It is currently unclear 
whether this cluster, which has a photometric age of 30 Myr, 
is associated with the knot, is merely a chance superposition, 
or is currently just passing near the knot. 

Whether or not this massive (sub) cluster is part of Knot 
S, it is clear that there is a significant age spread amongst 
the stellar populations present. This is not surprising, given 
the large length scales involved. In nearby galaxies, such as 
in the LMC, there is a clear correlation between the age 
difference between two star clusters and the physical dis- 
tance between them (Efremov & Elmegreen 1998). These 
authors interpret this as being due to the spatial extent of 
the progenitor giant molecular clouds with respect to the 
sound speed. Based on these authors' schematic representa- 
tion (made for the LMC), for a size of ~ 200 -500 pc, we 
would expect an age spread of 20-30 Myr, in good agree- 
ment with the observations. Based on this, we may expect 
W3 to have a similar, or even larger, age spread of stellar 
populations within it. 

Fellhauer & Kroupa (2005) have presented scaled N- 
body models of a complex of massive star clusters. In their 
models, the population of clusters is distributed in a cen- 
trally concentrated way; i.e. the clusters themselves follow a 
Plummer sphere with radius of ~ 100 pc. The stars within 
each cluster are in virial equilibrium, and the motions and 
distribution of clusters globally are also in virial equilib- 
rium. These authors show that (1) the clusters then merge 
in a relatively short timespan (a few Myr to a few 100 Myr 
based on simulations scaled to match W3) and (2) the re- 
sulting cluster features an extended envelope reaching out 
to ~ 1 kpc. 

While Knot S clearly has substructure within its halo, 
one of the predictions of the Fellhauer & Kroupa (2005) 
model is a dynamically cold core plus an extended profile, 
each of which can be described by a King or EFF profile. The 
resulting combined profile features a distinct bump where 
the two profiles join; i.e. the envelope has more stars than 
would be predicted from the extension of the inner profile 
into the outer parts. No such bump is seen in our profile 
of Knot S, although this cluster is probably too young to 
show such a behaviour, nor is such a bump seen in W3. The 
reason for this may be Fellhauer & Kroupa's (2005) assump- 
tion that all stars in the complex form in massive clusters 
which need to merge in order to disrupt. Possibly an initially 
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Figure 11. Comparison of the surface-brightness profile of Knot 
S in the Antennae galaxies and a YMC in NGC 6946 with that 
of W3. Knot S also has a power-law profile that extends to at 
least 400 pc, although its profile is shallower than that of W3. 
This difference could be intrinsic (i.e. present at formation) or 
alternatively due to evolution as W3 is ~ 400 Myr old, while 
Knot S is ~ 7 Myr old. A YMC in NGC 6946 is also shown 
(green) for comparison. This cluster also has an extended power- 
law envelope, as discussed in Larsen et al. (2001). Beyond the 
radius shown for this cluster, the profile becomes flat, due to the 
presence of a related stellar complex. 

smooth distribution of star formation with some clusters in 
it might dampen the resulting bump in the profile. Addi- 
tionally, the resulting profile might differ if the simulations 
began with a single dominant cluster at the centre of the 
complex surrounded by smaller clusters, instead of a number 
of massive clusters spread throughout the initial complex. 

6.2 Extended profiles in other YMCs 

While Knot S is the most massive known potential young 
counterpart to W3, being ~ 50 times less massive than W3 
was at birth (when taking only stellar evolution into ac- 
count), it is not the only extended cluster known. Larsen 
et al. (2001) studied the mass and profile of a ~ 15 Myr 
old YMC in the nearby face-on spiral galaxy NGC 6946. 
The profile of this cluster, shown as (green) diamonds in 
Fig. If II closely resembles that of W3 within the common ra- 
dial range. Using HST/WFPC2 images, Larsen et al. (2001) 
traced the profile out to 60 pc and found that it is well de- 
scribed by a power-law profile with a core radius of 1.2 pc 
and 7 = 2.1. 

The dynamical mass of this YMC is 1.7 x I0 6 M Q within 
65 pc of the cluster centre. Beyond this radius the profile 
becomes flat due to the presence of an underlying stellar 
association/complex that surrounds the cluster. It is pos- 
sible that after this surrounding association dissolves with 
time, the remnant profile will resemble that of W3 beyond 
65 pc. Like the younger Knot S, this YMC in NGC 6946 
also displays substructure in its envelope, with concentra- 
tions of stars and also smaller clusters. Larsen et al. (2002) 
find an age spread within the complex of ~ 30 Myr based 
on modelling the formation history of resolved stars. 

Closer-by, R136, the central star cluster in 30 Doradus, 
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is also known to feature an extended profile (e.g., Moffat 
et al. f994). Recently, Campbell et al. (20f0) have stud- 
ied Rf36 with the VLT using the Multi-Conjugate Adaptive 
Optics Demonstrator (MAD) to obtain images in the H and 
K passbands. The area surrounding Rf 36 is rich in gas and 
dust, and differential reddening effects can introduce spu- 
rious structure into the measured profile (e.g., Mackey & 
Gilmore 2003; Campbell et al. 20f0). Hence, imaging in the 
near-IR has the advantage of significantly reducing this bias. 
The profile of R136 from the Campbell et al. (20f0) study 
is shown as open (black) triangles in Fig. [TT] for radii below 
fO pc; it has a best-fit EFF-profile index of 7 = f .6. For radii 
larger than f pc we use the VXT/HAWK-I 7f-band profile 
derived from the integrated light, normalised to match the 
Campbell et al. profile at fO pc. 

This A'-band profile of Rf 36 matches the optical profiles 
of the YMC in NGC 6946, Knot S in the Antennae and W3 
in NGC 7252 out to ~50 pc. Moffat et al. (1994) report 
that the same power-law profile continues out to ~ f 00 pc 
(outside the HAWK-I field of view). Note that R136, along 
with its power-law envelope, makes up the dominant stellar 
population of 30 Doradus. 30 Doradus itself, over its full 100 
pc scale extent, is known to have an age spread within its 
stellar population of ~ f Myr (Walborn & Blades f997). 

In contrast, the Galactic cluster NGC 3603, which is 
similar to R136 in terms of radius, mass, density and age 
(Portegies Zwart, McMillan, & Gieles 20f0) appears to have 
a truncated profile, at least for the massive stars, at ~ f pc 
(Moffat et al. 1994). The youth of both NGC 3603 and Rf 36 
suggests that the differences between their profile types re- 
flect different initial conditions rather than different stellar- 
dynamical evolution (see Mai'z-Apellaniz 2001 for a discus- 
sion). 



6.3 A census of young massive stellar groupings 

Beyond these individual clusters, Mai'z-Apellaniz (2001) 
studied a sample of young (<20 Myr), relatively massive 
(> 3 x 10 4 Mq) stellar aggregates in nearby galaxies with 
HST images. He found that he could broadly split his sam- 
ple of 27 objects into three categories, with some overlap be- 
tween them: star clusters with "weak halos" (e.g., NGC 3603 
types), star clusters with "strong halos" (e.g., R136 types), 
and stellar aggregates that lack a distinct core (e.g., mas- 
sive OB associations). Within this categorization, R136, the 
YMC in NGC 6946, Knot S in the Antennae and W3, W30 
and W26 in NGC 7252 would all be classified as clusters 
with "strong halos" . ft is worth noting that all of the above 
galaxies with "strong halo" clusters also contain more com- 
pact "weak halo" clusters. 

It then appears that star clusters feature a range of 
outer profiles, with "strong halos" referring to clusters with 
shallow profiles, and "weak halos" to those with steep pro- 
files. 

As noted by Mai'z-Apellaniz (2001), the youth of many 
of the clusters and associations in his sample argues for the 
profile types to be due to the initial conditions during cluster 
formation rather than to dynamical processes after forma- 
tion. Do the profiles of clusters (extended, truncated, steep 
or shallow) vary systematically as a function of galactic loca- 
tion and environment? If so, what is the dominating factor? 



A systematic study of YMC profile types in nearby galaxies 
is required to answer these questions. 

6.4 Extended Envelopes in the NGC 7252 clusters 

While shallow and extended profiles have been found before 
for YMCs in a variety of galaxies (see above), the extended 
envelopes seen in the clusters of NGC 7252 are remarkable 
for at least two reasons. The first is their sheer extent. The 
envelopes of W26, W30 and W3 extend out to ~ f 50 pc, 
> 250 pc and > 500 pc, respectively. Only Knot S in the 
Antennae compares in spatial extent (see Fig. Iff [) . 

The second reason is the remarkable fact that these 
three clusters have been able to retain their extended en- 
velopes for 300 - 500 Myr. Knot S is only ~ 7 Myr old and 
hence may lose its extended envelope over the next several 
100 Myr. Similarly, one might expect the ~ 400 Myr old clus- 
ters in NGC 7252 to have lost their envelopes due to tidal 
truncation. That this did not occur for W3, W30 and W26 
is likely due to their extreme masses (> 5 x 1O 6 M0). In- 
terestingly, it is the lowest-mass cluster of these three, W26 
(~ 6 x fO 6 Mq) that shows possible signs of tidal truncation 
at ~150 pc (Fig. and 5 133)) . 



7 CONCLUSIONS 

We have presented an analysis of HST/WFC3 images of the 
galactic merger remnant NGC 7252. In particular, we have 
focussed on the structure and effective radii of the young 
massive clusters that formed during the merger. We have 
been able to resolve clusters down to ~ 2 pc, despite the 
64 Mpc distance to the host galaxy. The average effective 
radius of the resolved clusters in our sample is (Rea) ~ 6- 
7 pc, significantly larger than for YMCs in normal galax- 
ies, although we note that this may be caused by the fact 
that we cannot resolve the smallest clusters. We have found 
some evidence for a relation between the effective radius and 
mass for clusters with masses > 1 x 1O 6 M0, and also tenta- 
tive evidence for a relation between the effective radius and 
galactocentric distance, with clusters with larger projected 
distances having larger effective radii. 

For three of the brightest clusters, W3, W30 and W26, 
we have been able to measure surface-brightness profiles 
and constrain the structural parameters. We find that all 
three have power-law envelopes, well fit by EFF profiles to 
near the limit of measurability. These envelopes extend out 
to >500 pc, >250 pc, and ~ 150 pc, for W3, W30, and 
W26, respectively. Only the envelope of W26 shows a sign 
of being truncated (at ~ f50 pc), but whether this apparent 
truncation is real and tidal, or due to imperfect background 
subtraction, cannot presently be established with certainty. 
We note that it has been possible to trace the profiles out 
to these large distances mainly because of the brightness 
(high mass) of these three clusters and the relatively low 
and smooth background of the host galaxy. 

We have compared these three measured profiles to 
those of other YMCs found in the literature. A number of 
young clusters are known to also feature extended profiles, 
but only one — Knot S in the Antennae galaxies — extends as 
far out as those of the three NGC 7252 clusters. It is inter- 
esting to compare Knot S with clusters W3, W30 and W26, 



© 2012 RAS, MNRAS 000. fflfHl 



12 N. Bastian et al. 



due to the youth of its central region (~ 7 Myr). The enve- 
lope of Knot S contains a number of subclusters that may 
disrupt over the next few hundred Myr, leaving a smooth 
halo, as observed in the NGC 7252 clusters. The subclusters 
in Knot S show an age spread of ~ 30 Myr. Hence, if the 
progenitors of W3 and its siblings resemble Knot S, we may 
expect similar age spreads within them as well. 

The similarities between W3, Knot S, and other YMCs 
with extended envelopes suggest that the envelopes and their 
substructures are primordial in nature. Other very young 
clusters, such as NGC 3603 in the Galaxy, appear to form 
with a distinct truncation. We speculate that the differ- 
ence may stem from the gas distributions at the time when 
star formation begins: Gaseous cores distinct from their sur- 
roundings may form truncated star clusters, while core- 
cloud continua may form clusters with extended envelopes. 

In the process of our analysis, we have estimated the 
tidal radii of the three extended clusters in NGC 7252 from 
their estimated masses and projected galactocentric dis- 
tances, assuming that the clusters move on circular orbits. 
We find that these tidal radii are comparable to the ob- 
served extents of the profiles. This may explain why these 
very massive clusters (> 5x 10 6 Mq) were able to retain their 
extended envelopes for several hundred Myr. Other young 
clusters that display extended envelopes, such as R136 in the 
LMC or the YMC in NGC 6946, may lose their envelopes 
more quickly due to their lower masses and more disk-like 
environments. 
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